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SMALL  OBSTACLE  AND  APERTURE  THEORY 


I.  INTRODUCTION 

The  determination  of  the  scattering  parameters  (matrix)  or  equivalent  networks  for 
discontinuities  in  uniform  cylindrical  waveguide  usually  requires  the  evaluation  of  an  extremely 
complex  boundary  value  problem.  In  general,  solutions  are  not  easily  obtained  and  four  methods 
which  lead  to  approximate  solutions  are  classified  as 

(a)  The  variational  method 

(b)  The  integral  equation  method 

(c)  The  equivalent  static  method 

(d)  The  transform  method. 

\ 

The  reader  is  referred  to  the  Waveguide  Handbook,  Sec.  3.5,  for  a  thorough  discussion  of  the 
above  methods.  When  the  discontinuity  is  "small,”  however,  it  is  possible  to  deal  with  this  class 
of  problem  in  a  straightforward  and  simple  fashion.  A  small  discontinuity  is  defined  as  one 
which  is  located  "far"  from  the  waveguide  walls  and  whose  dimensions  are  small  compared  with 
wavelength.  When  these  conditions  are  satisfied,  the  scattering  of  the  electromagnetic  fields  by 
such  discontinuities  illuminated  by  a  known  excitation  is,  to  a  good  approximation,  independent 
of  the  effect  of  the  waveguide  walls  and  depends  only  on  the  geometry  of  the  discontinuity  and  the 
nature  of  the  illumination.  That  quantity,  which  is  solely  dependent  on  discontinuity  geometry 
and  excitation  field,  is  the  polarizability.  Since  the  discontinuity  is  small  compared  with  wave¬ 
length,  the  phase  of  the  excitation  field  is  essentially  constant  over  the  entire  discontinuity.  The 
polarizability  may  therefore  be  determined  from  static  field  considerations,  which  is  a  great 
simplification. 

It  is  the  purpose  of  small  obstacle  or  aperture  theory  to  combine  the  discontinuity  polariz¬ 
ability  with  object  position  in  the  waveguide  structure  and  the  incident  mode  fields  in  such  a  way 
that  the  scattering  matrix  and  equivalent  circuit  parameters  may  be  easily  determined.  It  should 
be  noted  that  not  only  does  this  approach  provide  a  simple  solution  to  otherwise  complex  prob¬ 
lems,  but  also  that  these  solutions  are  useful  in  regions  where  the  approximation  would  not  be 
expected  to  hold. 

The  results  which  are  presented  below  are  for  arbitrary,  lossless,  symmetric,  and  asym¬ 
metric  discontinuities  in  dominant  mode  uniform  cylindrical  waveguide. 
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II.  MODAL  ANALYSIS  AND  THE  REPRESENTATION  OF  FIELDS 

PRODUCED  BY  CURRENTS  AND  DISCONTINUITIES  IN  WAVEGUIDES 

The  basis  for  this  section  can  be  found  in  the  excellent  paper  by  Marcuvitz  and  Schwinger.^ 
Our  purpose  is  to  present  those  results  which  will  lead  directly  to  the  small  obstacle  or  aperture 
equations.  Some  typical  discontinuities  are  shown  below  in  Fig.  l(a-d). 
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(a)  (b)  (c)  (d) 

Fig.  1.  Typical  discontinuities  in  waveguides:  (a)  arbitrary  obstacle;  (b)  transverse 
aperture  or  obstacle;  (c)  longitudinal  aperture;  (d)  longitudinal  obstacle. 


It  is  a  fundamental  theorem  in  electromagnetic  field  theory  that  a  unique  solution  for  the 
electric  and  magnetic  fields  within  any  region  exists  provided  the  tangential  component  of  either 
the  electric  field  E  or  magnetic  field  H  is  specified  at  every  point  on  the  surface  bounding  that 
region  (including  discontinuity  points).  Thus,  in  the  steady  state,  the  electromagnetic  fields 
within  a  source-free  region  are  solutions  of  the  homogeneous  Maxwell  field  equations 


V  X  E  =  — jcjfjtH 

(la) 

V  X  11  =  jcoe  E 

(lb) 

where  an  e^wt  time  dependence  is  understood,  and  \i  and  e  are  the  permeability  and  dielectric 
constant  of  the  medium. 

For  an  obstacle  discontinuity  such  as  that  shown  in  Fig.  1(a),  the  region  of  interest  is  bounded 
by  the  waveguide  walls,  the  obstacle  surface,  and  the  far  "terminal"  surfaces.  A  far  terminal 
surface  (shown  by  a  dashed  line)  is  located  at  a  point  sufficiently  far  from  the  discontinuity  such 
that  only  dominant  mode  fields  are  considered.  The  higher-order  modes  generated,  being  eva¬ 
nescent,  may  be  considered  negligible  at  that  surface.  It  then  remains  to  solve  Eqs.  (1)  for  the 
electric  field  in  the  given  region  subject  to  the  boundary  values  of 


n  x  E 

on  the  guide  walls 

n  X  E 

on  the  terminal  surfaces 

H  x  n 

on  the  obstacle  surface 

(2) 

(n  =  outward  unit  normal). 

For  the  aperture  discontinuities  shown  in  Fig.  1,  there  are  two  regions  of  interest:  first, 
the  region  bounded  by  the  waveguide  walls,  the  aperture  surface  (shown  by  a  dotted  line),  and  the 
far  terminal  surfaces;  second,  the  corresponding  region  on  the  opposite  side  of  the  aperture. 

One  now  solves  Eqs.  (1)  for  the  magnetic  field  in  each  region  subject  to  the  boundary  values  of 

n  X  E  on  the  guide  walls 

H  x  n  on  the  terminal  surfaces 

n  x  E  on  the  aperture  surface  .  (3) 
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In  general,  it  is  extremely  difficult  to  solve  the  boundary  value  problems  for  the  regions 
illustrated  in  Fig.  1.  The  difficulties  arise  because  the  geometrical  shapes  of  these  regions  are 
perturbed  by  the  presence  of  the  obstacle  or  aperture  discontinuities  resulting  in  nonscparable 
coordinate  geometries.  Consequently,  it  becomes  most  convenient  to  invoke  the  concept  of 
"equivalence"  and  replace  the  discontinuity  structures  by  equivalent  electric  and  magnetic  sur¬ 
face  current  distributions.  This  approach  leads  one  to  consider  the  solution  of  the  inhomogeneous 
Maxwell  equations 


V  x  E  =  — jcjpH  -  M 

(4a) 

V  x  H  =  jweE  +  J 

(4b) 

in  a  discontinuity-free  region,  subject  to  elementary  boundary  conditions  where  the  fields  are 
produced  by  the  electric  and  magnetic  current  densities  J  and  M.  This  region,  a  perfectly  con¬ 
ducting  cylindrical  waveguide  of  arbitrary  but  uniform  cross  section,  is  shown  below  in  Fig.  2(a-b). 


Fig.  2.  Discontinuity-free  region:  (a)  transverse  view;  (b)  longitudinal  view. 

A  complete  representation  of  the  fields  for  such  regions  is  required  so  that  solutions  of  the 
inhomogeneous  field  equations  [Eqs.  (4)]. may  be  obtained.  Separation  of  the  total  field  solutions 
for  E  and  H  into  transverse  E^  and  Ht,  and  longitudinal  Ez  and  Hz  components  leads  to  a  great 
simplification.  The  inhomogeneous  field  equations  yield  for  transverse  fields 

d  E 

-  [11  +  VtVt]  •  (Ht  X  zo)  +  <Mt  X  zo)  (5a) 


3H 


sr  =  ni  +W  •  <zoxEt)  +  (zoXJt> 


(5b) 


and  for  the  longitudinal  fields 

z 

E  =  [V.  •  (H.  X  z  J  —  J  ] 

z  1  t  t  o  zJ 

z 

H  =  [V,  •  (z  X  E.)  -M  ] 

Z  JCJJJL  1  t  '  O  t  ZJ 

where  li  is  the  unit  dyadic,  Vt  =  V  -*zQ(a/az)  is  the  transverse  gradient  operator,  and 

M.  =  M.  -  —  V.  X  z  J 
t  t  JCJC  t  o  z 

J  =  J  -  -A-  V+  X  z  M 
t  t  JCJp  t  o  z 


(6a) 

(6b) 

(7a) 

(7b) 
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The  transverse  vector  field  equations  [Eqs.  (5)]  yield  the  scalar  transmission  line  equations 
if  all  transverse  vector  quantities  are  represented  in  terms  of  a  complete  set  of  orthonormal 
vector  mode  functions  (eigenfunctions).  The  representation  is  in  the  form  of  a  linear  superposi¬ 
tion  of  the  complete  set  of  eigenfunctions. 


and 


Et<r)  =  Z  v;(z)  e;<p> +  £  vr(z)  e;,(p) 

i  i 

Ht(r)  =  Yj  Jjtz)  h|(P)  +  £  Ij'tz)  h]'(p) 
i  i 


Mt(r)  =  Y  vj<z>  h-(p)  +  Z  vi'(z)  hi'(p) 

i  i 

jt(r>  =  y  qc**  +  Z  ii"<z>  ei'(p) 

i  i 


(8a) 


(8b) 


(9a) 


(9b) 


where  p  is  the  cross-sectional  vector  coordinate,  z  is  the  longitudinal  (symmetry)  coordinate, 
i  is,  in  general,  a  double  index,  and  ’  denotes  E-mode  and  "  denotes  H-mode.  The  longitudi¬ 
nal  field  components  are  written  as 


ju>eEz(r)  =  Y  E(z)Vt  •  h?  x  zq 
i 


(10a) 


jw(iHz(r)  =  Y  VV(z)  V  t  •  zQ  Xe!1 
i 


(10b) 


A  few  of  the  more  important  properties  of  the  eigenfunctions  (see  Refs.  1  and  2  for  more  complete 
descriptions)  are: 

E-Modes The  eigenfunction  e!(p)  may  be  derived  from  a  scalar  function 
. (p )  in  the  following  manner 

vt<pj(p) 

e.'(p)  =  —  V  uia) 

l  Kt. 

where  <p.(p)  is  a  solution  of 

(Vt2  +  k'2)  cp.(p)  =  0  (lib) 

subject  to  the  boundary  condition  <p.(p)  =  0  on  1  for  kt'.  ^  0;  [8cp.(p)/ 81]  =  0 

on  1  for  k'.  =  0.  In  addition, 
ti 

hj(p)  =  zq  X  e.r(p)  (11c) 

and 


e'  .(p)  =  —  j 
zi'K  J 


Y!k.,2<p.(p) 
1  tl^l'^ 


U€ 


(lid) 
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where 


Z!  =  1/Y!  =  <cj/w«  ;  «!2  =  k2-k’2  . 

II-Modes The  eigenfunction  h!'(p)  may  be  derived  from  a  scalar  func¬ 
tion  ip^(p)  in  the  following  manner 


VJAp) 

hi’(p)  = 

1  “ti 


where  ^.(p)  is  a  solution  of 


<v2+k”2)*i0»  =  ° 

subject  to  the  boundary  condition  [9^.(p)/an]  =  0  on  1.  In  addition, 

e! '(p )  =  h! '(p )  X  z 
1  K  l  o 


and 


h”<p)  =  -j 


Z."k?.,2i) 
1  tl  ^ 


:(P) 


C op 


where 


Z."  =  l/Y."  =  up/*."  ;  K."2  =  k2 -k;.'2 

l  '  l  r/  l  *  l  ti 

The  orthogonality  properties  of  the  vector  and  scalar  mode  functions  are 


(12a) 


(12b) 


(12c) 


( 1 2d) 


JjV  e*da 

s 

jl^jda= 

s 


6../k'2 
i y  ti 


hf  da  =  6.. 

3 

•  jlvi*-  vk;i2 

s 


(13a) 


(13b) 


where  <5...  is  the  Kronecker  delta  and  the  asterisk  denotes  complex  conjugate. 

Substitution  of  Eqs.  (8)  into  Eqs.  (5)  yields  the  following  generic  form  for  z-dependence  of  the 
fields  for  E-  or  H-modes: 


dV.(z) 

—  ,  =  j*r.Z.I.(z)  +  v.(z) 

dz  J  i  i  r  i 

(14a) 

dl.(z) 

- -3 —  =  i«-.Y.V.(z)  +  i.(z) 

dz  J  l  l  iv  i' 

(14b) 

Since  the  vector  mode  functions  form  a  complete  orthonormal  set, 
i.(z)  may  be  found  by  a  simple  inversion: 

the  source  terms  v^(z)  and 

Vi(z)  =  Mt<r)  *  (P)  da 

(15a) 

S 
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(15b) 


ij(z)  =  Jt(r)  •  e?  (p)  da 


S 

Substitution  of  Eqs.  (7)  into  Eqs.  (15)  obtains,  after  some  vector  algebra, 

h*  (p)  da 


/.(z)  =  Mt(r)  •  hr  (P)  da  + 


rz„  x  v*Jz(r> 


,ot“  ;■ 

a  s 

rn  pprzxv.M(r)i 

ii(z)  =  33  Jt(r)  • e* (p)  da  +  33  1  °  z  j  •  et (p)  da 


(16a) 


(16b) 


The  expressions  V.J  and  V  M  in  the  integrand  imply  that  J  and  M  are  differentiable  functions. 
I  Z  l  z  z  z 

Since  there  are  times  when  these  functions  are  not  differentiable  (e.g.,  point  sources),  the  form 
of  the  expression  is  changed  by  having  the  operate  on  the  mode  function  which  is  always  well 
behaved,  utilizing  Green's  theorem, 


if  v«  -  ,,a'  da  =  <£  f(A  •  n)  dl  -  ^  Vtf  •  Ada  +  ff[v.  •  Ada  (17) 

S  CSS 

where  1  is  the  bounding  contour,  n  is  the  unit  normal  to  1,  J ^  -  0  on  1,  h.  *  n  =  0  on  1,  and 
noting  that 


Zo  x  VtJz(r)  •  hi  (p)  =  ei  (p)  •  vtjz(r) 

Zo  X  VtMz(r)  •  e*  (p)  =  -h*  (p)  •  V^r) 


the  following  expressions  for  v.(z)  and  i.(z)  result: 


rr  rr  J  (r)  v+  *  ef  (p) 

hflPlda-j]  j‘. 


‘i|zl  '  ff  Jtlrl  '  e>lda-ff 


M 


,<r)  Vt  •  h*(p) 


da 


da 


(18a) 

(18b) 

(19a) 

(19b) 


S  S 

Since  V  *  e.(p)  =  jcoeZ.e  .(p)  and  V  •  h.(p)  =  jupY.h  .(p)  from  Eqs.  (11)  and  (12),  Eqs.  (19)  yield 
l  1  1Z1  XX  x  z  x 


v.(z)  =  M(r)  •  h*  (p)  da  +  Z*  J(r)  •  ej.(p)  da 
S  S 

i.(z)  =  ^  J(r )  •  e*  (p)  da  +  Yf  JJ  M(r)  •  h*.(p)  da  . 

S  S 

The  expressions  for  the  modal  source  functions  v.  and  i.  are  still  general. 


(20a) 


(20b) 


III.  THE  SMALL  OBSTACLE  APPROXIMATION 

Application  is  now  made  to  the  situation  where  discontinuity  dimensions  are  small  compared 
with  wavelength.  In  this  situation,  the  currents  induced  on  the  obstacle  are  approximately  in 
phase  and  the  predominant  effect  can  be  related  to  static  charge  distributions  which  vary  accord¬ 
ing  to  e^+  (the  time  dependence  of  the  excitation).  In  first-order  small  obstacle  (or  aperture) 
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theory,  these  induced  currents  depend  only  on  the  discontinuity  geometry  and  the  excitation;  the 
effect  of  the  higher-order  multipoles  and  of  the  surrounding  waveguide  walls  is  neglected.  One 
obtains  for  the  induced  electric  and  magnetic  dipoles,  Qe  and  Qm,  respectively, 

Qe(r)  =  ej>  •  E^nc(r)  (21a) 

Qm(r)  =  pHl  •  lCnc(r)  (21b) 

where  E^nc(r)  and  lCnc(r)  are  the  incident  electric  and  magnetic  fields  evaluated  at  the  obstacle 
center;  f  and  1\i  are  the  static  electric  and  magnetic  polarizability  dyadics,  respectively.  In 
addition,  if  6(r)  [=  6(p)  6(z)],  the  Dirac  delta  function,  defines  the  physical  extent  of  the  induced 
dipoles,  the  equivalent  electric  and  magnetic  current  dipole  elements  are  given  by 

J(r)  =  J<5(r)  =  jwQe(r)  <5(r)  (22a) 

M(r)  =  M6(r )  =  ju>Qm(r)  <5(r)  .  (22b) 

Substitution  of  Eqs.  (22)  into  Eqs.  (20)  yields 

v.(z)  =  v.S(z)  =  h*  •  M<5(z)  +  Z*  J  ■  e*.  S(z)  (23a) 

l'  1  lO  1  zio 


i.(z)  =  i.6(z)  =  e*o  •  J6(z)  +  Y*  M  ■  h*ioS(z)  (23b) 

where  the  subscript  o  denotes  evaluation  of  the  eigenfunction  at  the  coordinates  of  the  obstacle 
center.  The  transmission  line  equations  [Eqs.  (14)]  may  then  be  written  as 

dV.(z) 

-  =  jK.Z.I.(z)  +  v.«(z)  (24a) 

dl.(z) 

-  -gj-  =  J^YjVjtz)  +  E6(z)  (24b) 

and  yield,  upon  integration  between  z  =  0  and  z  =  0+, 

-[V.«T)  -  Vi(0+)j  =  v.  (25a) 

-[L(0-)  —  Ii(0+)]  =  i.  .  (25b) 

The  above  equations  imply  the  modal  equivalent  network  shown  in  Fig.  3. 

ijlcr)  i1(o+)  1 3  -46~ltfl6Fl 


Fig.  3.  Modal  equivalent  network. 
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The  voltages  and  currents  (V.  and  I.)  on  the  transmission  line  shown  in  Fig.  3  are  clearly 
given  by 


^J*;7 


V  .(z±)  =  -  4  [±v.  +  Z.i. ]  e 
1  2  1  i  i  r 


,  1 

Mz*)  -  zr  [±i.  +  Y.v.j  e 
r  2  1  l  i  iJ 


(26a) 

(26b) 


The  quantities  V.(z*)  and  Ltz*)  are  the  scattered  voltage  and  current  amplitudes;  the  mode  gen¬ 
erators  v.  and  i.  are  related  to  the  incident  voltage  and  current  amplitudes.  The  voltage  and 

11  4 

current  at  any  point  on  the  transmission  line  is  given  in  terms  of  traveling  waves. 


-IK.Z  „  1K.Z 

me  J  l  .  ref  J  i 
V.(z)  =  V.  e  +  V.  e 

i  i  i 


(27a) 


inc  ~JKiz  ref  ^ Kiz 
Z.I.(z)=V  e  1  -V  e  1 

i  iv  l  i 


(27b) 


The  normalized  incident  and  reflected  waves  as  defined  by  the  scattering  matrix,  S  (b  cSa),  are 
related  to  the  voltage  and  current  in  the  following  manner: 


1 

a.  =  — 
i  2 


V.(z) 


Ii<z) 


(28a) 


bi  =  I 


V.(z) 


(28b) 


One  therefore  obtains  the  following  for  a^  and  b.  for  a  two-port  structure  if  the  incident  energy 
at  ports  1  and  2  are  considered  separately: 


V.  -1K.Z 
10  ^  J  1 

a . .  =  -  e 

li  — 


(29a) 


b,.  = 

li 


.  l/C.Z 

-  [—v.  +  Z.i.]  e  1 

2  1  11 


(29b) 


a_ .  = 
2i 


V  i*.z 
10  J  1 
-  e 


(29c) 


_ .  -j/c.z 

b2i  =  - —  [v;  +  Z.i.)  e  1 


^  ^ 


(29d) 


The  electromagnetic  fields  may  be  concisely  described  in  terms  of  three-dimensional  electric 

th 

and  magnetic  vector  mode  fields  associated  with  the  i  mode,  having  unit  voltage  and  current 
amplitude,  respectively.  This  leads  to  a  particularly  compact  form  for  the  scattering  param¬ 
eters.  These  vector  mode  fields  are 


&f*hr)  =  (e.(p)  ±  ezi(p)]  e**'* 


(30a) 
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(30b) 


X  } 
1 


(*), 


(r) 


[h.(p)  ±  hzi(p)]  e 


*J*iz 


where  e.(p),  ezi(p)>  h^p),  and  h^.fp)  have  been  defined  previously  and  the  unit  current  is  given 
in  terms  of  the  unit  voltage  as  — Y.  for  waves  traveling  in  the  — z  direction  and  +Y.  for  waves 
traveling  in  the  +z  direction.  The  electric  and  magnetic  dipoles  J  and  M  are  then  given  by  (the 
(±)  superscript  distinguishes  between  an  equivalent  dipole  excited  by  a  wave  from  the  right  or 
left,  respectively) 


.J(±)  =  jwe?  •  (31a) 

M(±)  =  ijwjjiY.ln  •  Kj***  .  (31b) 


Thus,  the  scattering  parameters  describing  the  behavior  of  a  small  obstacle  (or  longitudinal  aper¬ 
ture  in  the  sidewall)  of  a  two-port  structure  are  (note  that  =  Z*  for  propagating  modes,  and 
£.(±)*  _  g.fc).  ft.W*  _  K.^  for  transversely  bounded,  lossless  waveguide) 
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(32a) 

(32b) 

(32c) 

(32d) 


The  above  formulas  will  yield  the  scattering  parameters  for  small  obstacles,  longitudinal, 
transverse  or  of  arbitrary  shape  in  dominant  mode  cylindrical  waveguide.  The  scattering  param¬ 
eters  for  small  longitudinal  apertures  may  also  be  obtained  if  one  replaces  the  object  obstacle 

A  A 

polarizabilities  9  and  ft,  by  its  dual-aperture  polarizabilities  9  and  ft,  which  can  be  shown  to 
have  values  equal  to  one-quarter  the  magnetic  and  electric  obstacle  polarizabilities.  Thus, 


A 

9  = 


(33a) 


(33b) 


On  the  other  hand,  the  scattering  parameters  for  small  transverse  apertures  can  be  deter¬ 
mined  most  easily  from  the  equivalent  circuit  of  a  small  transverse  obstacle  of  the  same  shape, 
size,  and  location  via  Babinet's  principle.  This  approach  will  be  carried  out  in  Sec.  IV.  The 
equivalent  circuits  for  small  obstacles  will  now  be  deduced  from  the  scattering  parameters  de¬ 
rived  above. 

Any  two-port  (two-terminal-pair)  network,  such  as  our  small  obstacle  (or  aperture)  con¬ 
figurations,  may  be  represented  most  generally  by  either  a  Pi  (admittance  formulation)  or  a  Tee 
(impedance  formulation)  equivalent  circuit  [see  Fig.4(a-b)]. 

The  relationship  between  the  scattering  matrix  S  and  the  impedance  and  admittance  matri¬ 
ces,  2  and  *2),  for  a  circuit  is  expressed  below 

Z  =  %f_4  =  (l  +  S)(1  —  S)"1  .  (34) 
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(a)  (b) 


Fig.  4.  Equivalent  circuits:  (a)  Pi,  and  (b)  Tee. 
Thus,  for  the  Pi  and  Tee  networks  shown  in  Fig.  4(a-b), 


where 
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where 


A  —  1  S  .  .  ®o')l'S..S__  s,0 

z  11  22  11  22  12 


(36d) 


In  the  small  obstacle  (or  aperture)  limit,  it  is  easily  shown  for  symmetrical  (y^  y^; 


z'  .  =  z'  )  discontinuities  that 
11  22 
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For  small  obstacles,  it  follows  directly  from  Eqs.  (32),  (35),  and  (36)  dial. 


y; 


12  sM  t  (1  -s12) 


'!2  (l-s12)-Sll 


Upon  substituting  Eqs.  (30),  (32a),  and  (32c)  into  Eqs.  (38),  one  obtains 
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(38a) 

(38b) 


(39a) 

(39b) 


IV.  THE  SMALL  APERTURE  APPROXIMATION 


The  scattering  and  equivalent  network  parameters  for  small  apertures  in  an  infinitely  thin, 
perfectly  conducting  screen  transverse  to  the  symmetry  (longitudinal)  axis  may  be  obtained  di¬ 
rectly  from  the  small  obstacle  results  via  Babinet’s  principle,  the  bisection  theorem,  and  the 
duality  properties  of  the  Maxwell  field  equations.  It  should  be  first  realized  that  a  transverse 
obstacle  may  be  represented  by  a  shunt  equivalent  circuit.  This  fact  is  easily  demonstrated  by 
an  examination  of  Eqs.  (39).  Since  the  transverse  components  of  the  magnetic  polarizability 
dyadic  and  the  normal  component  of  the  electric  polarizability  dyadic  are  zero  for  an  infinitely 
thin  obstacle,  it  follows  that 


Z11  z12 


2y 
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(40a) 


Y'  =  -=rr —  =  2(y'  -y'  )  =  ~r~  =  JCJ  [(lY.h*.  •  »  •  h  .  +  tZ.e?  •  9  •  e.  ] 

obs  Z'  'Jll  J12  z'  J  i  zio  zio  i  10  ioJ 

obs  12 

Figure  5(a-b)  shows  the  equivalent  circuit  of  a  thin  transverse  obstacle. 


1 3  -46  -  U862 ] 


(40b) 


Fig.  5.  Equivalent  circuit  of  thin  transverse  obstacle: 
(a)  admittance  description;  (b)  impedance  description. 


It  can  be  shown5  that  the  impedance  and  admittance  for  the  shunt  obstacle  and  aperture  are 
related  in  the  following  manner: 


=  _L  v  ’ 

4  obs 


(41) 
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Therefore, 


Z!=  j  ~  [tZ.e*  •  %  •  e  .  +  pY.h*  •  9  •  h.  1 

ap  J  4  i  zio  zio  i  10  ioJ 


(42) 


or,  in  terms  of  the  aperture  polarizabilities  %  and  9  [see  Eqs.  (33)], 


Z'  =  jo.  [|i.Y.h?  •  )li  •  h.  t  tZ.e* . 
ap  J  i  io  10  i  zio 
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(43) 


The  scattering  parameters  S11  (  =  Is^  |e  j  and  s12  (  =  ls12 le  j 
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11  1  +  2Z 1 

ap 

and,  for  lossless,  reciprocal  discontinuities, 

2  i  1 2  2  1 1 

si2  =  -C  ~  lsn  I2)2  e 


J^12\ 

e  ]  are  then  defined  by 


(44a) 


(44b) 


All  equations  derived  herein  are  consistent  and  exact  in  the  limit  of  discontinuities  of  zero 

1  6 

extent.  However,  it  has  been  shown  *  that  these  equations  provide  useful  data  even  for  discon¬ 
tinuities  whose  dimensions  are  significant  fractions  of  a  wavelength. 


V.  SELECTED  EXAMPLES 

The  following  examples  are  given  to  illustrate  the  ease  with  which  the  properties  of  other¬ 
wise  difficult  boundary  value  problems  may  be  obtained  using  small  obstacle  (aperture)  theory. 


Example  1 


Consider  an  elliptic  aperture  in  rectangular  waveguide  with  an  incident  mode  (see  Fig.  6). 


Fig.  6.  Elliptic  aperture  in  rectangular  waveguide 
with  incident  mode. 

From  Eq.  (43), 
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Example  2 

Consider  a  narrow  centered  strip  in  rectangular  waveguide  operating  in  the  dominant 
mode  only  (see  Fig.  7). 
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Fig.  7.  Narrow  centered  strip  in  rectangular  waveguide 
operating  in  dominant  H1Q  mode  only. 


From  Eq.  (40b), 
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Therefore, 


Hence,  upon  integrating  from  x  =  0  to  x  =  a,  one  obtains 
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